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We study AB n miktoarm star block copolymers in the strong segregation limit, focussing on the 
role that the AB interface plays in determining the phase behavior. We develop an extension of 
the kinked-path approach which allows us to explore the energetic dependence on interfacial shape. 
We consider a one-parameter family of interfaces to study the columnar to lamellar transition in 
asymmetric stars. We compare with recent experimental results. We discuss the stability of the A15 
lattice of sphere-like micelles in the context of interfacial energy minimization. We corroborate our 
theory by implementing a numerically exact self-consistent field theory to probe the phase diagram 
and the shape of the AB interface. 



I. INTRODUCTION 



Not only are block copolymers promising materials for 
nano-patterned structures 0, Q{, drug delivery and 
photonic applications Q , but they are also the ideal sys- 
tem for studying the influence of molecule architecture 
on macromolecular self-assembly J5J • Because of the on- 
going interest in novel macromolecular organization, the- 
oretical predictions based on heuristic characterization of 
molecular architecture offer crucial guidance to synthetic, 
experimental, and theoretical studies. Though the stan- 
dard diblock copolymer phase diagram Q was explained 
nearly a quarter of a century ago, the prediction and 
control of phase boundaries is fraught with subtle phys- 
ical effects: weak segregation theory provides an under- 
standing of the order-disorder transition Q , strong seg- 
regation theory (SST) predicts most of the ordered mor- 
phologies [8j, and numerically exact, self-consistent field 
theory (SCFT) @ can resolve the small energetic differ- 
ences between a variety of competing complex phases. 

In previous work, we argued that in diblock systems, 
that as the volume fraction of the inner block grows, AB 
interfaces are deformed into the shape of the Voronoi 
polyhcdra of micelle lattice, and therefore, the free- 
energy of micelle phases can be computed simply by 
studying properties of these polyhedra. In particular, we 
predicted that as volume fraction of inner micelle domain 
grows the A15 lattice of spheres should minimize the free 
energy as long as the hexagonal columnar phase (Hex) 
did not intervene |lflj . We corroborated this prediction 
by implementing a spectral SCFT Q for branched di- 
block copolymers: in this paper we probe the regime of 
validity of our analytic analysis through both strong seg- 
regation theory and SCFT. Though there is extremely 
small variation in the energy between different interfa- 
cial geometries, so too is the variation in energy between 
different stable phases. Thus, we compare these two ap- 
proaches not only by the phase diagram but also through 
the details of the ordering in the mesophases. Since our 
original ansatz hinged on the (minimal) area of the in- 
terface between the incompatible blocks, we will focus 
strongly on the shape and structure of this interface. We 
will explore in detail the relationship between molecule 
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FIG. 1: The miktoarm star architecture for an AB4 diblock 
copolymer is shown in (a). Shown in the Pm3n unit cell, 
(b), are the AB interfaces for the A15 phase, extracted from 
SCFT results for n = 5 at = 40 and / = 0.349 (along 
the Hex-A15 phase boundary). The A15 unit cell contains 8 
micelles: 2 at the center and corners (as with BCC); and 6 
along the faces of the unit cube. 



architecture and the polyhedral distortion of the AB in- 
terface induced by the lattice packing of micelles to study 
hexagonal columnar phases. Our results motivate the 
search for a stable A15 phase which we find in SCFT. 

In order to render the sphere-like phases stable in 
comparison to the Hex phase, we are obliged to con- 
sider asymmetric diblocks; while symmetric, linear di- 
block copolymers with an A and a B block have equiva- 
lent "inside-out" bulk morphologies when the A volume 
fraction / is replaced with the B volume fraction (1 — /), 
copolymers with branched or otherwise asymmetric ar- 
chitectures have no such symmetry, and therefore, tend 
to favor morphologies with one of the two components 
on the "outside" a of curved AB interface (i.e. on the 
outside of micelles). Indeed, our previous SCFT analy- 
sis of branched diblocks is consistent with these findings. 
Because of the challenge of accounting for all the com- 
peting effects, in Section II we implement a full SCFT for 
diblocks with the AB„ architecture to explore the mean 
field behavior of mitkoarm melts. In Section III we de- 
velop a strong-segregation theory approach for the hexag- 
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onal columnar phase which allows us to parameterize a 
large class of configurations and to explicitly assess the 
accuracy of the unit-cell approximation (UCA), which as- 
sumes the lattice Voronoi cell to be perfectly cylindrical 
(or spherical for three-dimensional lattices). Our calcu- 
lation builds on the "kinked-path" calculation of Milner 
and Olmsted and allows us to explore the influence 
of the hexagonal micelle lattice on the cylindrical mor- 
phology. We find that the shape of the Voronoi cell of 
the lattice strongly influences the shape of the AB inter- 
face. In Section IV we compare the predictions of the 
full SCFT calculation to the SST calculation in order to 
assess the accuracy of the latter. In addition, we demon- 
strate how the SST results of Section III can be used to 
compute an accurate phase boundary for transitions be- 
tween lamellar (Lam) to Hex configurations. We briefly 
discuss the inverse phases (where the B blocks are on the 
inside) in Section V. Finally, we conclude in Section VI. 



II. SELF-CONSISTENT FIELD THEORY OF 
AB„ MIKTOARM STAR COPOLYMERS 

Approximate self-consistent field theory calculations 
have explored the mean field phase behavior of linear 
diblocks with asymmetric monomer sizes 1 121 which were 
confirmed through numerically exact SCFT |l3Hl4| . Mil- 
ner developed SST, applicable in the xiV — > oo limit (x 
is the Flory-Huggins parameter for A and B monomers 
and N is the degree of polymerization of the copolymers) , 
for melts of A„B m miktoarm star copolymers which also 
incorporates asymmetric monomer sizes [ill 1 1 5| . Gener- 
ally, the results of all of these calculations show that equi- 
librium morphologies which have blocks with stronger ef- 
fective spring constants (i.e. more arms or smaller statis- 
tical segment lengths) on the outside of curved interfaces 
are favored over a much larger region of the phase space 
than in the symmetric diblock case. The details of the 
calculation implemented here will be reported elsewhere 
as a specific case of more general SCFT calculation for 
multiply-branched diblocks 0|. The method is an ex- 
tension of Matsen and Schick's spectral SCFT calcula- 
tion for melts of linear 0, ^| and starblock copolymers 
[l7j . Given the space group of a copolymer configuration, 
the mean field free-energy can be computed to arbitrary 
accuracy. The results of these SCFT calculations are ac- 
curate to the extent that mean field theory is correct and 
composition fluctuations can be ignored. The contribu- 
tions of these fluctuations tends to zero in the N — * oo 
limit |3 , and therefore we can expect to capture the equi- 
librium results observed in the PS-PI miktoarm star ex- 
periments, for which n ~ 1,000 m miu mm. 

We consider a melt of fixed volume and number of 
copolymers, with each molecule composed of N total 
monomers. The volume fraction of the A-type monomer 
is /. In general, we could accomodate monomer asym- 
metry by allowing for two different statistical segment 
lengths, a a and as, for the A and B species, respec- 



tively. In this case, each statistical segment length can 
be scaled appropriately so that the physical "packing 
length", £a,b = P~abi ^ a AB> * s nxe d to the proper value 
for each of the two chemical species Thus, without 
loss of generality we define a common segment density 
for the two monomer types, po- Moreover, for our SCFT 
calculations in this section we will restrict ourselves to 
the case a a — &B- The asymmetry of the copolymers 
we study arises entirely from their architecture. Each 
copolymer is composed of one block of pure A monomer 
joined to n blocks of pure B monomer at a common junc- 
tion point (see Figure ^ (a)). 



A. Low to Intermediate Segregation 

We computed the full phase behavior for %A < 40 for 
n < 5. To achieve a numerical accuracy of 0.005% for 
our free-energy calculations we employ up to 712 basis 
functions. This allows a precision in the phase boundary 
calculations which is better that ±0.001 for / and ±0.01 
for X-N- The computed SCFT phase diagrams are shown 
in Figures and 

The first notable feature of these phase diagrams is 
that they are not symmetric about / = 0.5, as is the case 
for symmetric diblocks - phase boundaries are shifted to 
the right. This indicates that phases with the n B blocks 
on the outside of curved interfaces are favored. By adopt- 
ing an interface with non-zero mean-curvature, the con- 
figuration relaxes the blocks on the outside of the AB 
interface at the expense of an increase in stretching of 
the inner blocks |23j . Due to the additional asymmetry 
introduced by the molecular architecture, the stability of 
phases with the B blocks on the outside of highly curved 
interfaces is enhanced. Moreover, this effect is generally 
amplified by further increasing the number of B blocks in 
the molecule. If we look, for instance, at the boundary 
between flat and curved interfaces which separates the 
gyroid (Gyr) from the lamellar (Lam) phase, we see that 
for %A = 25 the transition occurs at / = 0.358, 0.495, 
0.518, 0.523 and 0.525 for n = 1, 2, 3, 4, and 5, respec- 
tively. Likewise, the "inverted" morphologies, with B 
blocks on the inside of curved interfaces, are suppressed 
because curving the interface inward towards the B do- 
main introduces excess stretching in these blocks. This 
pushes the phase boundaries to greater / for the inverse 
structures as well. 

We also note the appearance of a stable A15 phase 
of spherical micelles (see Figure (b)). The A15 phase 
is also observed in experiments and simulations of an- 
other soft molecular system, namely neat mixtures of low 
molecular weight dendrons [M EH EE El The stabil- 
ity of this phase has been attributed to the fact that 
the area of the A15 Voronoi cell is the minimal among 
the Voronoi cells of three dimensional periodic structures 
(for a given number density) [5E EH ■ We have argued 
that this minimal area makes the A15 lattice stable in 
the phase diagram of a 3-generation, multiply-branched 
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FIG. 2: Phase diagrams for n = 2 and n — 3. Dis labels re- 
gions where the melt is disordered. Stable regions of ordered 
phases are labeled: (Lam) lamellar; (Gyr) gyroid, Ia3d sym- 
metry; (Hex) hexagonal-columnar, pQmm symmetry; (A15) 
sphere phase, Pm3n symmetry; (BCC) body-center cubic lat- 
tice of spheres, Im3m symmetry; and (FCC) face-centered 
cubic lattice of spheres, Fm3rn symmetry |4q| . The circle 
marks the mean field critical point through which the system 
can transition from the disordered state to the Lam phase via 
a continuous, second-order phase transition. All other phase 
transitions are first-order. 

copolymer m the (extreme) limit that the AB in- 

terface adopts the flat shape of the Voronoi cell of the 
micelle lattice, the A15 phase is stable over other sphere 
phases, such as BCC and FCC. Thus, the appearance of 
the A15 phase in these phase diagrams is evidence of a 
highly distorted AB interface when B blocks compose the 
outer domains of micelles. Due to melt incompressibilty, 



FIG. 3: Phase diagrams for n = 4 and n = 5. Stable phases 
and critical points are labeled as in Figure [5] 

copolymer configurations must fill the lattice Voronoi cell 
so that the monomer density is uniform [30l | . The tension 
of the outer block chains stretching towards the corners 
of the cell pulls on the interface, distorting it towards the 
polyhedral shape of the Voronoi cell. Of course, the inter- 
facial tension of the AB interface and the tension of the 
inner blocks will frustrate this distortion, and the inter- 
face will actually distort into a rounded polyhedra of the 
same symmetry as the Voronoi cell. In the next sections 
we will return to this point and will show, nonetheless, 
that the approximation of flat-faced AB interfaces im- 
proves as the number of outer arms grows. 

Indeed, from Figures [21 and [3] it is clear that increasing 
the asymmetry increases the stability of the A15 phase; 
the window of A15 stability increases further in both the 
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larger / and smaller / directions. Since increasing n, the 
number of B blocks per molecule, stabilizes sphere-like 
phases over cylindrical phases for larger / the A15-Hex 
boundary should move to larger /. Moreover, as n grows 
the effective spring constant in the outer blocks grows as 
n 2 : for fixed number of B monomers (1 — f)N, there are 
n entropic springs with spring constant [a 2 (l — f)N/n]~ 1 . 
This increased resistance to stretching will lead to a 
greater polyhedral distortion of the AB interface. This 
interfacial distortion is apparent in the oblate shape of 
the A domains of A15 phase in Figuren(b), correspond- 
ing to the oblate shape of Voronoi cells of lattice sites on 
the faces of the Pm3n unit cube. For larger values of n 
interfaces approach the shape of Voronoi cell at smaller 
values of / and thus the BCC-A15 transition occurs at 
lower /. Of course, for the inverted micelles we expect the 
opposite to be true; the tension of the highly stretched 
inner domains will prefer a spherical interface, ignoring 
the shape of the Voronoi cell. Therefore, we would not 
expect that the inverse A15 phase is stable at large / and 
we do not find it within SCFT. 

At low segregations, near the order-disorder transi- 
tion (ODT), the predicted phase behavior is dramati- 
cally altered from the case of linear diblocks. Generally, 
the ODT is shifted up to higher ^JV, or lower tempera- 
ture. For symmetric diblocks the critical point, indicat- 
ing a mean field, second-order disordered (Dis) to lamel- 
lar transition, occurs at \N — 10.495 and / = 0.5 0. 
The critical points of the miktoarm star copolymer melts 
are shifted, for example, %A = 13.47 and / = 0.559 for 
n = 2 and \N = 13.94 and / = 0.588 for n = 3. The up- 
ward X-/V shift in the ODT is consistent with other SCFT 
calculations for conformationally asymmetric copolymers 
[l2l IT3L Il4| and indicates that chain fluctuations near the 
ODT are systematically altered by molecular asymmetry. 
Finally, we note that neither the cubic double-diamond 
nor the hexagonal-perforated lamellar phases are stable 
in these systems as with symmetric, linear diblocks. 



B. Strong Segregation 

In order to compute the strong segregation phase 
behavior of AB„ miktoarm star copolymer melts, we 
compute the Dis-BCC, BCC-Hex, and Hex-Lam phase 
boundaries at \N ~ 100 |3l|. Due to the numerical 
difficulties of considering the Gyr and A15 phases for 
large %V we do not include these phases in our stabil- 
ity analysis. Additionally, we ignore the window of sta- 
ble closed-packed spheres, FCC, which is predicted by 
SCFT calculations to occur near the ODT. Since we ex- 
pect that the free-energy differences between the sphere 
phases will be relatively small at these strong segrega- 
tions 0, |n| our calculation should capture the change 
from spheres to cylinders. Likhtman and Semenov devel- 
oped a general SST calculation to assess the stability of 
bicontinuous phases in the limit \N — > 00. Their calcu- 
lation shows that the Gyr phase is unstable when chain 
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FIG. 4: SCFT results for AB n miktoarm stars at large \ N - 
For our calculation, the asymmetry parameter e = hb = n, 
the number of B blocks per molecule. Triangles, A, depict 
Dis-BCC transitions; open circles, o, depict BCC-Hex tran- 
sitions; and diamonds, o, indicate Hex-Lam transitions. All 
boundaries are computed at \N — 100 with the exception of 
the low-/ BCC-Hex and Hex-Lam boundaries for n — 3, 4 and 
5. For n = 3 these boundaries are computed at \N = 80, and 
for n = 4 and 5 these boundaries are computed at \N — 60. 
Dark circles, a, indicate equilibrium results from experiments 
on PI-PS melts [%| . S, C, and L label spherical, cylindrical, 
and lamellar morphologies, respectively. 



fluctuations are ignored 32]. It is not known whether 
the bicontinuous Gyr phase is stable for finite but large 
XN even within the well studied SCFT phase diagram of 
symmetric, linear diblocks |33l |. Recent experiments on 
melts of fluorinated PI-PEE diblocks indeed suggest that 
this morphology is an equilibrium phase for \N ^100 
|34j . and we expect that the same may be true for more 
asymmetric copolymers, such as these miktoarm stars. 
Although we cannot prove its stability with our SCFT 
calculations, we expect that the Gyr phase appears at 
compositions intermediate to the stable Hex and Lam 
regions. 

In Figure 0] we compare our SCFT results for AB m 
copolymers to the experimental results on polyisoprene 
(PI)-pob/stvrene (PS) A 2 B 2 AB 2 19], AB 3 Hf, and 
AB5 |2l], [23] miktoarm star copolymer melts. Milncr 
showed that SST phase boundaries depend only on /, 
and the asymmetry parameter, e = (ns /nA){a,A/ &b)i 
where ub and ua are the number of A and B chains 
per molecule |15| . Since we are considering symmetric 
monomers, a a = as, "a = 1, and ub — n for our AB n 
copolymers, we have simply e = n. 

We see from Figure 0] that SCFT very accurately 
captures the observed phase behavior of miktoarm star 
copolymer melts. Therefore, ignoring composition fluc- 
tuations, as SCFT dictates, does not alter the phase be- 
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havior at this level of segregation. Moreover, as the data 
suggest, the effect of the molecular asymmetry seems to 
saturate for n > 3, and phase boundaries do not change 
significantly as a function of / for further asymmetry 
n. It has been suggested that the saturation of phase 
boundaries occurs when the spheres or cylinders of the 
inner block reach close packing, i.e. / ss 0.90 for Hex 
and / w 0.68 for the body-centered cubic (BCC) lattice 
|2g . In the next section we employ SST to explore how 
the symmetry of the micelle lattice frustrates their self- 
assembly, and moreover, how this frustration is related 
to the saturation of phase boundaries for e> 1. 



III. STRONG-SEGREGATION THEORY: 
"KINKED-PATH" REVISITED 

The appearance of a stable A15 phase suggests that 
the shape of the AB interface is strongly affected by the 
lattice symmetry for highly asymmetric diblock melts. 
Tension in the outer block chains maintains a uniform 
outer domain thickness and consequently distorts the in- 
terface into the shape of the Voronoi cell. However, pre- 
vious studies have demonstrated that for symmetric AB 
diblocks the shape of the interface, and therefore, the 
calculation of the free-energy are insensitive to the shape 
of the Voronoi cell HH|< The free energy is dom- 
inated by the tension of the AB interface, and so the 
minimal-area cylinder is favored. Because of this, a unit- 
cell approximation (UCA) is often taken for the domain 
shape of the micelles. For instance, in the columnar and 
spherical phases, the Voronoi cell is approximated by a 
perfect cylinder or sphere, respectively. However, pack- 
ing these cylinders or spheres into a space-filling lattice 
leaves voids in the interstices which is incompatible with 
the incompressible melt state. Therefore, the UCA can 
only provide an estimate for the free-energy of these mor- 
phologies. In fact, this estimate is a lower-bound to the 
true free-energy since distorting the round, approximate 
unit cells can only raise the free energy either by stretch- 
ing the outer block or distorting the interface. As our 
results in the last section suggest, the "packing frustra- 
tion" 35] between the surface tension and the stretching 
is highly dependent on molecular asymmetry. It is also 
unlikely that the UCA can capture the saturation as a 
function of the number of blocks, n; the close-packing 
limitation cannot be captured in the UCA since the in- 
ner volumes can fill 100% of the approximate unit-cell 
without over lap [13. 

In order to quantitatively explore the role of packing 
frustration in miktoarm star copolymer melts, we build 
upon the SST calculation of the Hex phase free energy 
for miktoarm star copolymers Olmsted and Mil- 

ner developed a "kinked-path" ansatz for the extension 
of the copolymer chains, which allows the configuration 
to maintain a cylindrical AB interface while satisfying 
the constraints of melt incompressibility (see Figure EJ. 
Here, we extend that calculation allowing the interface 




(a) (b) 

FIG. 5: In (a), a schematic representation of the kinked-path 
ansatz for a hexagonal unit cell with a circular interface di- 
viding A and B domains, depicted with a dashed line. Chains 
in the domain extend radially from the micelle center, while 
chains in the B region stretch away from the wall to the cor- 
ners of the cell, modifying the shape of the outer domain 
wedges. Symmetry dictates that the B blocks extend radially 
along the six mirror planes of the unit cell, shown as arrows. 
We need only consider the bending of chain paths within the 
right-triangular wedge shown in (b). Both angles, 6 and 0, 
are measured with respect to the y-axis. The A blocks extend 
along Ra(#), while the B blocks extend along D, and R(O) 
is the vector from the center of the micelle to the edge of the 
outer domain wedge at the wall of the unit cell. 



to adopt a more general class of interfaces, allowing us to 
systematically explore the effect of packing frustration in 
this morphology. 

Within SST we ignore chain fluctuations, in addition to 
the composition fluctuations absent from the full SCFT. 
We assume that each chain extends only along its clas- 
sical trajectory, and these paths must be consistent with 
the constant monomer concentration of the melt state. 
Following jllj we divide the hexagonal unit cell into in- 
finitesimal wedges which extend along the direction of 
the chain paths. Due to incompressibility we must have 
wedges with a volume fraction of A domain / and a vol- 
ume fraction of B domain (1 — /). Unless the AB in- 
terface adopts the same shape as the Voronoi cell, the 
chains, and consequently the wedges must bend towards 
the cell corners in order to distribute the volume of the 
chains evenly. 

Each chain starts with an A block that extends from 
the center of the micelle radially to the AB interface at 
an angle, 9 (see Figure (b)). This path extends along 
the vector, R^^), which parameterizes the interface as 
function of 9. From that point on the interface, the outer 
blocks extend to some other point, at angle O, on outer 
wall of the unit cell, along the vector D(0). The vec- 
tor which extends from the center of the micelle to the 
outer wall of the unit cell at O is given by R(0). We 
can parameterize the kinking of the wedges by the func- 
tion 0(9) which maps the chain position at the inter- 
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face to the position at the wall of the cell. Therefore, 
the extension of the B portion of the wedge is given by 
D(0) = R(0(0))-R A (0). 

The hexgonal unit cell is parameterized for < < 
7r/6 by 



R(9) = i? c 



2\/3 



secOr(O) , 



(1) 



where ttRq is the cross-sectional area of the micelle and 
r(/3) = sin/3x + cos/3y is the radial unit vector. Be- 
cause the configuration is strongly segregated and the 
inner domain is composed only of A monomers, we con- 
sider a class of interfaces which encloses a constant area 
fraction, /, of the total micelle area. In particular, we 
would like to consider a class of interfaces which vary 
from area-minimizing and uniformly curved (circle) to 
stretch-minimizing and polygonal (hexagon) . Such a one 
parameter family of interfaces is: 



R A (9;a) = f 1/2 R J{l-a)+a[ 



( 7r 



sec 2 9 r(0) (2) 



When a = this is the circle and when a = 1 it traces 
out an affinely shrunken version of the hexagonal unit 
cell JJ|. For all other values of a E [0, 1], the shape in- 
terpolates the two extremes (Figure EJ. By construction 
the enclosed area of the interface is JttRq for all a. 

We compute the free energy as a function of a and 
minimize in search of a variational ground state. We 
expect the ground state to favor the minimal area, round 
interface (a = 0) in the low / limit because differences 
in stretching to the walls and corners of the Voronoi cell 
are small on the scale of the radius of gyration of the 

1 /2 

outer blocks, aUl — f)N/n\ . In the other limit, / — > 
1, differences in chains stretching towards the corners 
and walls becomes large on this length scale, and the 
hexagonal interface (a = 1) will be preferred. We refer to 
the former case as the round-interface limit (RIL) and the 
latter as the polyhedral- interface limit (PIL). Of course, 
there is no guarantee that the true ground state belongs 
to this class of configurations. Nonetheless, minimizing 
the free energy over a will provide a variational upper- 
bound on the ground state free energy which is lower 
than either that from the straight-path calculation or the 
circular interface calculation. 

The essence of the kinked-path approximation is that 
each wedge must locally satisfy the volume constraint for 
the diblock so that: 



/ [R(0) - R A (9)} x — [R(0) + R A (9)} = 

(1 - /)R A (0) x Ar a (0) ( 3 ) 

Given the parameterizations introduced in Eqs. Q and 



(0) we find that 
tan (6(6»)) = 



(l-«)* + a(^)taa*-^(&^)sii 



-2V3) 



tt ( R A (a,6 
V3\ Ro 



cos 9 



(4) 



where R A (a,9) — \Ra(oi,9)\. It is easy to verify that 
tan (0(0)) has the appropriate limits: when a = 0, Eq. 
(0} reduces to the results of ref. |n|, and when a = 1, 
tan (0(6*) J = tan#, which is the straight-path ansatz. 

The free energy has three parts, arising from the in- 
tcrfacial energy, the stretching of the A blocks and the 
stretching of the B blocks: 



F = Fi, 



F A + F E 



(5) 



The interfacial energy per molecule (in units of ksT) 
is given by the effective surface tension, 7 ~ x 1 ^ 2 H3j 
times the area of the AB interface divided by the number 
of chains per micelle: 



2NJ 1 / 2 1 A 



Po 



Ra 



(6) 



where ttRq is the cross sectional area of the micelle and, 



A(a) = 




(7) 



is the area of the AB interface measured in units of 
the area of circular interface enclosing the same volume. 
We compute the stretching free-energy usin g th e self- 
consistent parabolic brush potential for melts [38|, yij as 
is shown in ref. Illl. For the inner A domain we have 



F A = 



lQNa\ 



I(a)R 2 , 



where 



12 rJ 6 d9R\{9, a) 
l{a) = - 



2nR*p 



(l-a 2 ) + 



5tt 



(8) 



(9) 



is the reduced stretching moment of the A domain, intro- 
duced in ref. The stretching term due to the outer 
B domain is given by 



F B 



3ir 2 n 2 B 



8(1 - /) 2 JVa; 



-S(aJ)R 2 



(10) 



where the stretching moment of the B blocks is given by 
the integral 



S(a, /) = 
d9T>{9)- 



N A (0) 



dR A 



(19 



+3N(0(0)) 



ds 



D 2 {9) 
~^Rj 



(11) 
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FIG. 6: Plots of the optimal interfacial shape parameter, a, 
versus composition, /, for e = 1, 2, 3, 4, and 5. On the right 
are the interface shapes given by Eq. for a = 0, 1/3, 2/3, 
and 1, shown as dashed contours. Inset on the upper left is 
the hexagonal unit-cell showing the AB interface (dashed line) 
extracted from the SCFT calculation for miktoarm stars with 
n = 3, \N = 80 anQl / — 0.54. The minimum and maximum 
distances from the center of the micelle to the interface, r m i n 
and r m ax respectively, are labeled. These are used to compute 
an estimate of the shape parameter for the SCFT calculations, 
a', defined in Eg, lib! . These values of a' are depicted as 
diamonds, <>, for n = 2 and \N = 100, and triangles, A, for 
n = 3 and = 80. 



with 



and 



N(e(0)) 



N A (0) 



dR 



(W 



Ro 



2VS 



dtan.8 



(10 



= R A {0)r 



(10 



dR A 
— -( 

dJS 



y , (12) 



(13) 



with 8 — cos #x — sin 6y. 

Combining these terms we have the total free-energy 
per molecule and minimizing over the dimension of the 
micelle, Ro, we compute the micelle free-energy per chain 
in units of fc^T: 



F Hex (a) =F A(a) 



2/3 



2fl(a) 12efS(aJ) 
e (l-/) 2 



1/3 



(14) 

where e = (jib /nA){o-Al ' a s) and Fq is proportional to 
(%iV) 1 / 3 and is independent of composition pT| . 

The inner domain stretching, 1(a), and the interfa- 
cial area, A(a), are minimized by circular interface at 
a = 0, while the outer domain stretching, S(a,f), is 
always minimized for a > 1 (though we are only con- 
sidering a € [0,1]). As / — > 1 the only finite energy 
configuration is the hexagonal interface with a = 1. Min- 
imizing over alpha for 0.15 >/ > 0.65 and 1 > e > 5 



we find the optimum energy interface configurations and 
show the results in Figure In the low / limit a round 
interface is favored due to the relatively lower tension in 
the B chains. For larger /, the stretching term of the B 
chains l|14l) begins to dominate, distorting the interface 
towards the hexagon. Because of the 0(e 2 ) increase in 
the effective spring of the outer chains, it follows that for 
larger asymmetry the onset of this polyhedral distortion 
occurs at lower compositions /. 

For high molecular asymmetries, e > 3, and A vol- 
ume fractions larger than / ss 0.65, we find equilibrium 
shapes with a > 1. In this case, the interface is bowed 
in, away from the walls, and pulled out towards the cor- 
ners. However, because it is the only finite energy con- 
figuration, each solution approaches a — 1 as / — > 1. 
We expect that the true ground state prefers an inter- 
face which asymptotically approaches a hexagonal shape, 
rather than "overshooting" a = 1 before returning to 
the / — ► 1 limit. The bowing inward of the interface 
is due to the oversimplified representation of the inter- 
face which does not allow the corners of the interface to 
relax smoothly towards the corners of the Voronoi cell. 
Nevertheless, this simple representation of the interface 
shape proves sufficient to capture the effect of the pack- 
ing frustration introduced by the hexagonal symmetry of 
the micelle lattice. 

In order to make contact with the results of SCFT and 
experiment, we need to estimate the value of a for an 
arbitrary interface. We introduce a numerical measure 
of the distortion from the circular shape. For copoly- 
mer melts in the large \N limit the ratio of the inter- 
facial thickness to the domain size scales as (x^) -2 / 3 
[23J; thus, in the strong segregation limit the interface 
becomes infinitely thin on the scale of the domain. We 
find the location of the interface by the contour for which 
the volume fraction of A is 0.5. One of these contours 
is shown in the inset of Figure EJ Using this contour we 
measure the minimum and maximum distance from the 
center of the micelle to the interface, r m i n and r max , re- 
spectively. As a measure of the distortion from a circular 
interface, we introduce 



(15) 



which is similar to the normalized amplitude of the 
cos(6#) modulation of the interface [3!|. Since 6 = 
for the circular interface and 5 = 7 — 4y3 m 0.0718 for a 
hexagonal interface, we may define 



(16) 



7-4V3 



In the two limits a = = a' and a = 1 = a' and it 
can be shown that {a' — a) /a is never greater than 5% 
for the class of interfaces parameterized by Ra{ch, 0). In 
particular, (a' - a) /a » 0.05,0.04,0.03,0.02,0.01, and 
for a = 0,0.2,0.4,0.6,0.8, and 1, respectively. 

By using SCFT for linear diblock copolymers, Matsen 
and Bates showed that a' w 0.006 for / as large as 0.34 
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IV. SELF-CONSISTENT FIELD VS. 
STRONG-SEGREGATION 




n QfS r i i i , i i i i i i i i 

0.2 0.3 0.4 0.5 0.6 

/ 

FIG. 7: Plots of the free-energy of the optimal configurations 
in units of the free-energy of the PIL, straight-path configu- 
ration, F H ex = F H ex(a min ) I F He x(a = 1). The dashed line 
shows the straight-path upper bound, where the Voronoi cell 
and interface are hexagonal. The dotted line shows the UCA 
lower bound, where the unit-cell and interface are circular. 



[35j , and they concluded that there is negligible distortion 
so that the effect of packing frustration can be ignored for 
the Hex phase. However, we find that even for the n = 1 
diblock case, interfacial distortion becomes appreciable 
f° r / ^ 0.5. While the Hex phase is unstable at these 
compositions for AB diblocks, experiments for melts of 
ABC triblocks yield stable Hex phases with a composite 
AB inner domain volume fraction of about 0.65 [4(j. For 
these triblocks, the interface separating the outermost 
C domain and the B domain is found to be very nearly 
hexagonal, demonstrating that packing-frustration in the 
Hex phase is indeed amplified as the volume fraction of 
the inner domain is increased. Indeed, this frustration 
can be relieved by adding C-type homopolymer which 
aggregates in the corners of the hexagonal unit cell, al- 
lowing the BC interface to relax to uniform curvature 

We plot the free-energy of the optimal configuration as 
a function of composition for the various asymmetries in 
Figured For lower compositions the free-energy is nearer 
the UCA lower bound. As the composition increases the 
free-energy approaches the PIL upper-bound. Clearly, 
for large e and large / the free-energy of our configu- 
ration is best approximated by the straight-path upper 
bound This suggests that for asymmetric copoly- 

mers, the AB interface of the Hex micelle is significantly 
deformed from the uniformly-curved configuration over 
compositions ranges where this phase competes for sta- 
bility. We review the consequences of this conclusion for 
the phase calculation in the next section. 



We compare the results of our SST calculation for 
asymmetric miktoarm star copolymer melts to the "more 
exact" results of our SCFT calculations for the Hex 
phase. Since the SST approximation results from an 
asymptotic expansion of the full SCFT partition function 
in the X N -> oo limit 0, we expect that SST and SCFT 
results should agree provided that we have used the ap- 
proximately correct shape of the AB interface. Strictly 
speaking we should consider our SST results to be an up- 
per bound on the free-energy of the Hex phase since we 
did not minimize over all possible ground state interface 
and chain configurations. Nevertheless, we argue that 
this bound is sufficiently close to the true ground state 
to capture the phase behavior of asymmetric miktoarm 
stars. 

We use the correspondence between a and a' for 
i?vi(a, 0) as a basis for direct comparison of the predicted 
values from SST and the "measured" values of a 1 from 
SCFT. A comparison is made in Figure for n = 2 and 
ri = 3. Despite our simplified representation of interface 
shape, Ra{ol 1 9) 1 and the crude means of measuring the 
equilibrium SCFT shape parameter, the predicted a and 
computed a 1 compare remarkably well. In particular, the 
SCFT results confirm that increasing / and e increases 
the equilibrium shape distortion. There is a systematic 
tendency in SST to overestimate the hexagonal distor- 
tion for low inner domain volume fractions. This can 
be attributed to the inadequacy of our interface param- 
eterization, Ra(cc,0), at low distortions where the AB 
interface is likely best represented as a constant radius 
with a superposed cos(6#) modulation j3f| . Nevertheless, 
the agreement at larger / can be taken as evidence that 
we have accurately accounted for packing frustration for 
our SST Hex phase calculation. 

The effect of interfacial distortion on the location of 
the predicted Lam-Hex phase boundary for miktoarm 
star copolymers is shown in Figure |S| There we show 
(with dashed lines) the calculated phase boundaries us- 
ing both the lower- and upper-bound approximations 
which assume that the interface and Voronoi cell are ei- 
ther both circular (i.e. the UCA) or both hexagonal (i.e. 
the PIL), respectively pj|. Note that there is no dif- 
ference between the round and polyhedral Voronoi cell 
for the Lam phase, and therefore, no approximation is 
necessary for this geometry. Since the UCA calculation 
underestimates the Hex free-energy it overestimates the 
composition at which the Hex phase becomes unstable 
against the Lam phase. Likewise, the upper-bound calcu- 
lation underestimates the location of the Lam-Hex phase 
boundary. We expect, therefore, that the true Lam-Hex 
SST phase boundary should lie between these two ap- 
proximate boundaries. 

Though the upper-bound of the round-interface limit 
is close to the UCA lower-bound for symmetric diblocks 
, the latter cannot be used to predict the phase behav- 
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FIG. 8: The dashed line shows the predicted Lam-Hex phase 
boundaries computed using the PIL, upper-bound (leftmost 
boundary) and UCA, lower-bound (rightmost boundary) ap- 
proximations for the SST Hex phase. The solid line is the 
Lam-Hex boundary computed using the SST results of Sec- 
tion III. The dotted lines are the Hex-BCC phase boundaries 
by comparing the PIL, upper-bound (leftmost) and UCA, 
lower-bound (rightmost) SST results for those morphologies. 
For comparision the \N ~ 100 SCFT results for the Lam- 
Hex boundary are shown as diamonds, o, and the BCC-Hex 
boundary are shown as open circles, o. 



ior for asymmetric diblock melts in SST, as experimental 
results suggest [22( ■ This discrepancy is illustrated in Fig- 
ure |H1 while SCFT predicts that the effect of increasing 
asymmetry saturates for e > 3, the UCA predicts that 
the phase boundaries move to larger f as e increases. 
The solid line in Figure OH shows the SST phase Lam-Hex 
phase boundary computed with the results of the previ- 
ous section. For e < 2 the AB interface is nearly undis- 
torted and circular; thus, the free-energy is best approx- 
imated by the lower-bound, UCA calculation. Hence, 
our Lam-Hex phase boundary follows the UCA phase 
boundary for low asymmetry. However, for e > 2 the ef- 
fect of interfacial distortion becomes important, and the 
free-energy is closely approximated by the upper-bound, 
PIL calculation. Therefore, for large asymmetry our SST 
Lam-Hex boundary tends toward the hexagonal-interface 
approximation phase boundary. This transition from the 
round to the hexagonal interface approximation is also 
reflected in the SCFT boundary. Indeed, our SST and 
SCFT Lam-Hex phase boundaries agree well in this en- 
tire parameter range (1 < e < 5). 

Similarly, we see that Hex-BCC boundary computed 
from SCFT tends to follow the UCA boundary for 
e < 2 and follows the polyhedral-interface approxima- 
tion boundary for e > 2, confirming our picture of the 
effect of packing frustration on asymmetric copolymers. 



Again, we know that interfacial distortion must be im- 
portant along the phase boundary between sphere and 
columnar phases because the A15, with its minimal in- 
terfacial area, is the favored sphere phase there. It would 
be interesting to pursue a one-parameter class of inter- 
faces that interpolated between a spherical and polyhe- 
dral AB interface. Such a calculation would not only pre- 
dict the spherical-cylindrical micelle phase boundary but 
it would also elucidate the transition from BCC spheres 
with nearly spherical interfaces at small / to A15 sphere- 
like micelles with polyhedral interfaces for larger /. 

We can also see that within the polyhedral-interface 
approximation calculations, the effect of increasing asym- 
metry beg ins to saturate for e > 5, as seen in exper- 
iment [22| . Again, this stems from the relative e 2 in- 
crease in stretching energy of the B blocks over the 
A blocks. In the large e limit the A block stretch- 
ing can be ignored, and l|14|) can be approximated by 

F Hex (a min ) « F [l2eM(l) 2 5(l,/)/(l-/) 2 ] 1/3 , where 
we have assumed that a m i n ~ 1< The same is true 
for the free-energy of other morphologies, namely Fx ~ 
^ e 1/3 ^(/), where X = Lam, A15, BCC, etc. and 
&x{f) is some function specific to the morphology which 
depends only on composition Because all free- 

energies scale the same way with e, the location of the 
transition between any two phases depends only on com- 
position and is independent of e in this limit. It should 
be noted that the effect of increasing asymmetry satu- 
rates for the UCA bounds as well, though this occurs 
very close to / = 1 for both the Lam-Hex and Hex-BCC 
boundaries. The location of phase boundaries computed 
using the polyhedral-interface approximation saturates 
near / w 0.50 and / « 0.75, for the Hex-BCC and Lam- 
Hex boundaries, respectively. 

While the agreement between, experiment, our SCFT 
calculations and our SST analysis for e < 5 is improved 
over previous efforts, the SST results appears still to over- 
estimate the effect of increasing asymmetry. We should 
note, however, that SST is an asymptotic expansion of 
the full SCFT near %iV — > oo, and at finite values of %iV 
we need to assess the importance of higher order correc- 
tions. It is known that lowest order, strong-segregation 
results are OKxN) 1 / 3 ] and that the leading corrections 
which can distinguish between phases are 0(1). One of 
these corrections is associated with fluctuations of chain 
ends and junction points and are thus proportional to 
(riB + tia — 1) J3- The other is associated with a prox- 
imal layer near the interface, where the stretching en- 
ergy of the chains deviates from the predictions of the 
parabolic brush potential |43j. At x^ ~ 100 this cor- 
rection can be as high as 20%, and will surely effect the 
relative free-energies of competing phases. Further errors 
arise from the implicit assumption of a parabolic chemi- 
cal potential used to compute the free-energy associated 
with molten polymer brushes [38l l39| . It has been shown 
that this approximation is equivalent to allowing for a 
negative chain end-density near the surface of a convex 
brush; nevertheless, it has been shown that relaxing the 
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constraint of a non-negative end-distribution hardly per- 
turbs the free-energy of a two-dimensional convex brush 
|44j . However small, this parabolic chemical-potential 
underestimates the constrained stretching free-energy of 
domains on the outside of cylindrical and spherical do- 
mains. A correction to the stretching energy of B do- 
mains at large asymmetries of 1% leads to a correction of 
the full free-energy of order 0.3%. While this seems small, 
the difference between upper- and lower-bound SST free- 
energies is only 3.6% and 6.8%, for the Hex and BCC 
phases, respectively. Since even small corrections to the 
free-energy lead to significant changes in the predicted 
phase behavior, we expect approximations of this order 
may be relevant. 



V. THE INVERSE PHASES 

We have applied the SST analysis here to the Hex 
phase where the A domain composes the interior of the 
micelle and the multiple B blocks compose the outer do- 
main. These asymmetric configurations have the "stiffer" 
blocks on the outside of micelle, hence, they are respon- 
sible for appreciable interfacial distortion. However, we 
could also use our SST results to determine the phase 
behavior of the inverse micelles, where the multiple B 
blocks composed the inner domain. Assuming that all 
B chains extend radially from the center of the micelle 
to the AB interface, we can compute the inverse micelle 
free-energy by the substitutions e — > e _1 and / — > (1 — /) 
in Eq. Q. 

We find, however, that such a calculation produces 
overestimates of the free-energy. Moreover, the assump- 
tion that the B chains stretch only radially is inconsistent 
with SCFT results for the inverse phases. Figure (a) 
shows a real-space A monomer distribution for the Hex 
phase with n = 4, x N = 100 and / = 0.3. Figure H (b) 
plots the distribution of A monomer along a horizontal 
line extending across the micelle for / = 0.3 and the in- 
verse phase at / = 0.835. Clearly, both configurations 
are strongly segregated with the volume fraction of the 
A monomer, <Pa(z), constant outside the inner domain 
(cj) A (z) = for / = 0.3 and (j> A {z) = 1 for / = 0.835). 
However, while the volume-fraction of B component is 
zero inside the / = 0.3 micelle, the volume fraction A 
component is never zero inside the inverse micelle at 
/ = 0.835. This shows that in the inverse micelle, the 
miktoarm star junction points are not strictly confined 
to a narrow AB interface. Rather, junction points are 
located somewhere within the inner domain and some 
of the chains must bend back towards the outer domain 
instead of stretching radially towards the center of the 
micelle. Thus the simple SST analysis of Section III 
does not apply to inverse phases since the assumption 
that junction points are confined to the AB interface is 
clearly violated. In order to account for this we would 
parameterize the distribution of junction points in the 
micelle, and minimize the free-energy over the distribu- 
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FIG. 9: Shown in (a), a real-space reconstruction of the A 
monomer volume fraction profile, <j>A (z), for the Hex phase 
at n = 4, xN — 100, and / = 0.3. High concentrations of 
A monomer are shown as bright red regions, and low concen- 
trations of A monomer (high concentrations of B) are shown 
as dark blue regions. In (b), the A monomer volume frac- 
tion profile along the horizontal line in (a) for xN ~ 100 
and / = 0.3, the solid line, and / = 0.835 (the inverse Hex 
phase), the dashed line. The equilibrium length scales, D, 
are l.MN 1/2 a and 1.40iV 1/2 a, for / = 0.3 and / = 0.835, 
respectively. 



tion. Moreover, we would need to find the optimal con- 
figuration of B blocks, bending inwards and outwards in 
the micelle. Such a calculation introduces many more 
degrees of freedom; minimizing over these is likely very 
difficult. Since SCFT provides numerically exact results 
for even very large segregations, it is a much better ap- 
proach to finding the inverse micelle configurations to 
determining the phase behavior. 
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VI. CONCLUSION 

We have implemented SST and SCFT calculations 
which elucidate the coupling between the molecular 
asymmetry of diblock copolymers and packing frustra- 
tion introduced by the micelle lattice. In particular, we 
find that the shape of the AB interface is highly sensitive 
to molecular asymmetry, in the case of our miktoarm 
star copolymers, the number of B blocks per molecule. 
For linear diblocks, or otherwise nearly symmetric ar- 
chitectures, we show that the AB interface of micelles 
in the Hex phase is rather insensitive to the shape of 
Voronoi cell and maintains a nearly constant mean curva- 
ture shape. However, for AB„ miktoarm star copolymers 
with n > 3 the interface is highly distorted towards the 
hexagonal shape of the Voronoi cell in regions where the 
Hex phase competes for stability. The effect of this dis- 
tortion is to shift the predicted Hex-Lam phase boundary 
to lower compositions than is predicted by approximating 
the Voronoi cell as a perfect cylinder. 

As a consequence of the importance of interfaces, we 
expect the A15 phase of spherical micelles to be stable for 
highly asymmetric copolymers due to the minimal area 
of its Voronoi cell amongst lattices in three-dimensions. 
This prediction is borne out by our numerical SCFT cal- 
culations which compute the full phase behavior of AB„ 
miktoarm star copolymers for n < 5. Indeed, we find 
that as molecular asymmetry is increased, the stability 
of the A15 phase is enhanced indicating that the polyhe- 
dral approximation for the AB interface is more valid as 
the copolymer architecture becomes more asymmetric. 

It is worth noting that the A15 phase of spherical mi- 



celles has yet to be identified experimentally in melts of 
highly asymmetric diblock copolymers. Pochan, Gido 
and coworkers report the appearance of an equilibrium 
cubic phase of spheres in PS-PI AB2 copolymer melts 
in regions of the e — f phase space where we might ex- 
pect A15 to be stable 19] . However, they note that 
the small-angle X-ray scattering data from these melts 
cannot distinguish between a simple-cubic or BCC mi- 
celle lattice; nor does the data rule out the possibility 
of the A15 lattice which contains all of the reflections 
of the BCC lattice. Experiments on sphere phases of 
more asymmetric miktoarm star copolymers in param- 
eter ranges where A15 should be stable tend to yield 
poorly ordered arrangements of micelles 0]. Because 
the A15 and BCC phase are nearly degenerate (with a 
free-energy difference on the order of 0.2%), the system 
may be kinetically trapped in some sort of glassy inter- 
mediate state in these parameter ranges. Despite these 
difficulties, however, both the robustness of the minimal 
Voronoi cell area principle and the SCFT calculations 
strongly indicate that the A15 phase of spherical micelles 
is an equilibrium phase of copolymer melts of sufficient 
molecular asymmetry. 
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